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Ground State of the Three-Band Hubbard Model 
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The ground state of the two-dimensional three-band Hubbard model in the oxide superconductors is 
investigated by using the variational Monte Carlo method. The Gutzwiller-projected BCS and SDW 
wave functions are employed in search for a possible ground state with respect to dependences on elec- 
tron density. Antiferromagnetic correlations are considerably strong near half-filling. It is shown that 
the d-wave state may exist away from half-filling for both the hole and electron doping cases. Overall 
structure of the phase diagram obtained by our calculations qualitatively agrees with experimental in- 
dications. The superconducting condensation energy is in reasonable agreement with the experimental 
value obtained from specific heat and critical magnetic field measurements for optimally doped samples. 
The inhomogeneous SDW state is also examined near 1/8 doping. Incommensurate magnetic structures 
become stable due to hole doping in the underdoped region, where the transfer t pp between oxygen 
orbitals plays an important role in determining a stable stripe structure. 

PACS numbers: 74.20.-z,74.25.Dw,71.10.Fd. 

I. Introduction bard model based on perturbative calculations such as 

generalized RPA treatments. 26-28 In perturbative treat- 
ments of the one-band and three-band Hubbard models, 
the spin fluctuations induced by the on-site Coulomb in- 
teraction promote anisotropic pairing correlations. QMC 
evaluations with some constraints due to the fermion sign 
problem are against a possibility of superconductivity in 
the three-band Hubbard model. 25 



In order to investigate the mechanism of supercon- 
ductivity (SC) in cuprate high-T c superconductors, 1 we 
examine the ground state of the two-dimensional three- 
band Hubbard model for CuC"2 planes which are con- 
tained usually in their crystal structures. It is believed 
that the Cu02 plane contains the essential features of 
high-T c cuprates. 2,3 It is not an easy task to clarify the 
ground state properties of the 2D three-band Hubbard 
model because of strong correlations among d and p elec- 
trons. We must treat the strong correlations properly to 
understand the phase diagram of the high-T c cuprates. 
The quantum variational Monte Carlo method (VMC) 
is a tool to investigate the overall structure of phase di- 
agram from weak to strong correlation regions. In this 
paper we investigate possible ground states in the three- 
band Hubbard model for Cu02 plane by employing VMC. 

Superconductivity in the one-band Hubbard model has 
been studied by numerical 4-13 and analytical 14-19 calcu- 
lations. The three-band Hubbard model has also been in- 
vestigated with intensive efforts recently. 20-29 The exact 
diagonalization computations for the three-band model in 
early stage of high-T c research supported a possibility of 
superconductivity by showing that holes can bind in small 



systems 



30,31 



It is also reported that the attractive inter- 



action works for both the <i-wave and extended-s wave 
channels based on finite temperature quantum Monte 
Carlo (QMC) simulations. 21 It has been shown recently 
that one can predict finite T c for the three-band Hub- 



In order to investigate the possibility and origin of su- 
perconductivity, the recent work by Kondo is important 
where it has been shown that the d-wave state has lower 
energy than the normal state for small U by employing 
the perturbation theory in U for the one-band Hubbard 
model. 32 This indicates that the ground state is super- 
conductive with ci-wave symmetry for small values of U. 
We can expect that this also holds for the three-band 
model. 33 It is then natural to expect that the d-wave 
state is stable for finite U unless there occurs some or- 
dering in the ground state. Among several possible long- 
range orderings, antiferromagnetic one should be exam- 
ined because the state with antiferromagnetic ordering is 
considerably stable near half-filling. In fact, according to 
VMC work for the one-band Hubbard model, the antifer- 
romagnetic (AF) energy gain is larger than the SC energy 
gain by almost two order of magnitude near half-filling. 
Then the competition between SC and AF states is very 
severe for the SC state. 12 ' 13 The SC region for the one- 
band Hubbard model is considerably restricted and a pos- 
sibility of pure superconducting state is very small. 12 A 
similar feature has been obtained by VMC evaluations for 
Ud = oo three-band Hubbard model 34 where antiferro- 



f 



magnetic region extends up to 50 percent doping and the 
d-wavc phase exists only in the infinitesimally small re- 
gion near the boundary of antifcrromagnctic phase. Thus 
VMC results performed recently are consistent with the 
constrained path QMC calculations 25 in the sense that 
a possibility of d-wave phase for the one-band Hubbard 
model and Ud = oo three-band Hubbard model is small 
at present, although an attractive interaction works for 
d-wave pairing. 

We expect that the antiferromagnetic region will shrink 
for the three-band Hubbard model if we adjust param- 
eters contained in the model. The parameters of the 
three-band Hubbard model are given by the Coulomb 
repulsion Ud, energy levels of p electrons e p and d elec- 
tron e<i, and transfer between p orbitals given by t pp . A 
purpose of this paper is to investigate the property of 
antiferromagnetic state and a competition between anti- 
ferromagnetism and superconductivity for finite Ud based 
on the three-band model following ansatz of Gutzwillcr- 
projected wave functions. 

It has also been argued that holes doped in the anti- 
ferromagnetically correlated spin systems induce incom- 
mensurate spin correlations in the ground state for the 
one-band Hubbard model 35-40 and three-band model 41 
within the mean field approximation. In the mean- 



H = 
+ 
+ 
+ 



where 

V = U d J2dj^dj l d il . (2) 

i 

x and y represent unit vectors along x and y directions, 
respectively. p\ ±& / 2 a an< ^ Pi±x/2,a denote the operators 

for the p electrons at site Ri ± x/2. Similarly p\±y/ 2 a 
and Pi±y/2,a are defined. Other notations are standard 
and energies are measured in units of t dp - For simplicity 
we neglect the Coulomb interaction among p electrons. 

We consider the normal state, BCS and SDW wave 
functions with the Gutzwiller projection. These types of 
functions are standard wave functions and well describe 
the ground-state properties with several long-range order- 
ings. They have been investigated intensively for the one- 
band Hubbard model. n ~ 13 ' 51 ~ 55 In Refs. 51 ' 55 it has been 
discussed that they can be improved systematically by 
operating correlation factors c~ XHa c~ aV . For the model 



field treatment the energy scales appear to be extremely 
large compared to values for real materials. Recent 
neutron-scattering experiments revealed incommensurate 
spin structures 42-49 developed at low temperatures and 
at low energies. The static incommensurate struc- 
ture was reported on LSCO samples: I^-zSr^CuCU, 
Lai.6-a;Ndo.4Cu04 and La2- x Sr x Ni04+j / . The incom- 
mensurate magnetic peaks have been also reported for 
YBa 2 Cu307_,5 by the inelastic neutron-scattering exper- 
iments. This type of inhomogeneous state may possibly 
provide a key concept to resolve the anomalous proper- 
ties of high-T c cuprates in the undcrdoped region. We 
will examine a possible phase of incommensurate states 
for the three-band Hubbard model by variational Monte 
Carlo method. 

The paper is organized as follows. In the next sec- 
tion the wave functions are presented. The SC state and 
uniform SDW state are discussed in Section III and a sta- 
bility of incommensurate state is examined in the subse- 
quent section. A summary is given in the last section. 

II. Hamiltonian and Wave Functions 

The Hamiltonian is given a s 25 ' 34 ' 50 



(1) 



shown above they are written as 

^u^Pg II a LlO), (3) 

\k\<k F ,rj 

V-sc = PqPn. Y[(u k + Wfe4 T aI fei )|0), (4) 

k 

^sdw = Pg II ^ll°>' ( 5 ) 

\k\<kp ,cr 

where a^ a is the linear combination of dk a , p x ka and 
Vyka constructed to express an operator for the lowest 
band of a non-interacting Hamiltonian in the hole picture. 
For t pp = 0, is expressed in terms of a variational 
parameter £ p — id'- 



*dpy^[ d L(P»+a/2,q J rPi+y/2,a ~ Pi-x/2,* ~ Pi-y/2,a) + h.C.] 
ia 
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H + V, 
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where w x k = 2td p sin(k x /2), w y k = 2td P sm(k y /2), w k 



w 



yk 



|V2 an d E k = [(e p - e d ) 2 /A + wl} 1 / 2 . The 



Fourier transforms of d- and p- electron operators are 
defined as 



Pxka 



Pyka 



1 



N l/2 



1 



_/yl/2 /_^ i 1 j i+x/2<j 



t 



Jk-(Ri+x/2 



Nl/2 



+y/2a 



Jk-(Ri+y/2} 



(7) 



(8) 



(9) 



where N is the total number of cells which consist of d, p x 
and p y orbitals. Coefficients Uk and Vk, appearing only 
as a ratio, are given by the BCS form: 



Uk 



z k + (e k + A 2 y/ 2 - 



(10) 



for = £fc — M where €k is the energy dispersion for the 
lowest band. Pq is the Gutzwiller projection operator 
for the Cu d site and P/v E is a projection operator which 
extracts only the states with a fixed total electron num- 
ber. The SC order parameter is assumed to have the 
following d x 2_ y 2- and extended s-wave form: 



d A fc = A s (cosfc x — cosky), 
S* Ak = A^cosfcr + cosky). 
Equation (4) is written as 



N c /2 



(11) 

(12) 
(13) 



The wave function given by eq.(13) agrees with 

4>BCS = PGY[(u k +v k al 1 J_ kL )\0), (14) 

k 

in the thermodynamic limit. For the commensurate 
SDW state (3ka is given by a linear combination of two 
wave numbers k and k + Q for the commensurate vector 
Q = (n, it). We can also investigate the incommensurate 
SDW state with incommensurate vector Q = (w ± 2ir8, ir) 
by diagonalizing the Hartree-Fock Hamiltonian with an- 
tiferromagnetic long-range order. The system sizes are 
given by 6 x 6 and 8x8 for the projected BCS wave func- 
tion and 16 x 4, 24 x 6, 32 x 8, 40 x 10 and 16 x 16 for the 
incommensurate SDW states. Our calculations are per- 
formed with the periodic and the antiperiodic boundary 
conditions for the x- and y-direction, respectively. This 
set of boundary conditions was chosen so that does 
not vanish for any k-points possibly occupied by elec- 
trons. 
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FIG. 1. Energy per site (E nor 



E)/N of the SDW state 



as a function of hole density S for t pp =0.4 and Ud = 8. From 
the top, e p — ed — 3, 2, 1.5 and 1. The results are for 6x6, 
8x8, 10 x 10 and 16 x 12 systems. Antiperiodic and peri- 
odic boundary conditions are imposed in x- and y-direction, 
respectively. Monte Carlo statistical errors are smaller than 
the size of symbols. Curves are guide for eyes. 
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FIG. 2. Energy per site {E normal - E)/N of the SDW state 
as a function of hole density S for t pp — 0.0, 0.2 and 0.4 where 
e p - e d = 2 and U d = 8. The results are for 6 x 6, 8 x 8, 10 x 10 
and 16 x 12 systems. Curves are guide for eyes. 



The expectation values are calculated following the 
standard Monte Carlo procedure by using the Metropolis 
algorithm. In the process of finding a minimum of en- 
ergy, we should optimize many parameters included in 
the wave functions. For such purpose we employ corre- 
lated measurements method to reduce the required cpu 
time. 56 

III. Condensation energy and phase diagram 

First, let us discuss the SDW phase near half- filling 
by evaluating the ground-state energy for optimized pa- 
rameters g, e~ p — £d and AF order parameter Aaf- We 
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FIG. 3. (a) Energy per site (E norma i - E)/N of the SDW 
state as a function of hole density 8 for Ud = 8, 12 and 20 
where e p — = 2 and t pp = 0.2. (b) Antiferromagnetic order 
parameter as a function of hole density for Ud = 8, 12 and 20 
where e p — ed = 2 and f pp =0.2. The results are for 6x6, 
8 x 8, 10 x 10 and 16 x 12 systems. Curves are guide for eyes. 



set e p = throughout this paper. It is expected that 
holes introduced by doping are responsible for the disap- 
pearance of long-range antiferromagnetic ordering. 57-59 
We show the SDW energy gain AE$dw in Fig.l as a 
function of doping ratio for several values of e p — ed- 
AEsdw increases and the SDW region becomes large 
as e p — Cd increases. The figure 2 shows the SDW energy 
gain for several values of t pp , where AEsdw is reduced 
as tpp increases. In Figs. 3(a) and 3(b) the dependence on 
Coulomb repulsion Ud is shown; the SDW phase extends 
up to 30 percent doping when Ud is large. Then it follows 
that the SDW region will be reduced if e p — Cd and Ud de- 
crease or tpp increases. In fact, Fig. 4 shows the boundary 
of SDW phase in the t pp -5 plane for Ud = 8 where 8 is 
the hole density and negative density indicates electron 
doping. Compared to the calculations for Ud = oo the 
SDW region is reduced greatly. 34 

Next, let us turn to the projected-BCS wave function, 
where the Gutzwiller parameter g, effective level differ- 
ence e p — id, chemical potential fx and superconducting 
order parameter A s are considered as variational param- 
eters. In Fig. 5 we show the energy as a function of A s 
where t pp — 0.0, Ud — 8 and e p — ed = 2 and doping ratio 
is given by 8 = 0.111 for (a) and 8 = 0.333 for (b). The 
d-w&ve superconductivity is most stable among various 



FIG. 4. Boundary of the SDW state in the plane of t pp -8 for 
e p — £d = 2 and 1. We set Ud = 8. 

possible states such as isotropic s-wave and anisotropic 
s-wave pairing states. The squares in Fig. 5 denote the 
values for the normal state, which are estimated indepen- 
dently by using an alternative Monte Carlo algorithm. 
The finite SC energy gain indicates that the attractive 
interaction works for <i-wave pairing. 

The SC energy gain (which is called the SC condensa- 
tion energy in this paper) is also dependent on e p — ed, 
as is shown in Fig. 6 for t pp — 0.2, Ud — 8 and 8 = 0.111 
on 6 x 6 lattice. This shows a tendency that the SC con- 
densation energy increases as e p — ed increases, which is 
consistent with calculations for Ud — oo. 34 It is noted 
that the dependence on e p — ed for the SC energy gain is 
rather weak compared to the SDW energy gain. We also 
note that the SC energy gains for Ud = 8 are mostly of 
the same order of those for Ud = oo. 34 

From the calculations for the SDW wave functions, we 
should set e p — ed and Ud small so that the SDW phase 
does not occupy a huge region near half-filling. In Fig. 7 
we show energy gains for both the SDW and SC states for 
Ud = 8, tpp = 0.2 and e p —ed = 2, where the negative delta 
indicates the electron-doping case. Solid symbols indicate 
the results for 8x8 and open symbols for 6 x 6. For 
this set of parameters the SDW region extends up to 20 
percent doping and the pure d-wave phase exists outside 
of the SDW phase. The d-wave phase may be possibly 
identified with superconducting phase in the overdoped 
region in the high-T c superconductors. 

The superconducting condensation energy obtained by 
our calculations is estimated as E con d — 0.0005^^ = 
0.75meV per site in the overdoped region near the bound- 
ary of SDW phase from the difference between the min- 
imum and the intercept of the E/N — A s curve with 
the vertical axis, where we set td p = 1.5eV as esti- 
mated from cluster calculations. 60-62 We have also es- 
timated E con d from several experiments such as specific 
heat or critical field measurements for optimally doped 
samples. They are given as 0.17 ~ 0.26meV from specific 
heat data 12,63,64 and 0.26meV from critical magnetic field 
value -ff 2 /87r. 12 ' 65 Our value is in reasonable agreement 
with the experimental data as was already shown for the 
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FIG. 5. Ground state energy per site as a function of A s 



on 6 x 6 lattice 
5 — 0.333 and t~ 



for (a) S = 0.111 and t p 



0.0, and (b) 



0.0. Parameters are given by Ud = 8 and 



£p — £d — 2 in units of tdp- Squares denote the energies for the 
normal state evaluated independently. 
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FIG. 6. Superconducting (circles) and antiferromagnetic 
(squares) energy gains per site as a function of e p — for 
t pp = 0.2 and Ud = 8 on 6 x 6 lattice. 



FIG. 7. Condensation energy per site as a function of hole 
density S for t vp = 0.2, e p — = 2 and Ud = 8. Circles, squares 
and diamonds denote the energy gain per site in reference 
to the normal state energy for d-wave, SDW and extended-s 
wave states, respectively. Solid symbols are for 8 x 8 and open 
symbols are for 6x6. Curves are guide for eyes. 



Hubbard model where the SC energy gain in the bulk 
limit is given by 0.00117t/site=0.59meV/site. 54 This 
agreement between the theoretical and experimental con- 
densation energy is highly remarkable. We expect that 
this value is not far from the correct value according to 
the evaluations for improved wave functions, 51 where it 
was shown that the energy gain is not changed so much 
due to multiplicative correlation factors e~ ATHo e~ ArV . 
We cannot estimate the SC condensation energy in the 
underdoped region because the SDW state is more stable 
than d- wave state and the SC condensation energy is not 
available experimentally due to a loss of entropy in the 
underdoped region. 63 

The phase structure obtained by our calculations 
agrees well with the available phase diagram indicated 
by experiments qualitatively, which means that a large 
SDW phase exists in the underdoped region and there 
is a G?-wave superconducting phase next to SDW phase 
in the overdoped region. Our calculations for electron- 
doping case predict d-wave symmetry away from half- 
filling, which is consistent with recent experiments on 
Ndi. 85 Ce .i5CuO4 



66 



IV. Incommensurate antiferromagnetism with spin 
modulation 

In this section let us discuss the underdoped region 
where the SDW state is significantly stable as shown in 
the previous section. Let us note that the SDW state 
can be possibly stabilized further if we take into ac- 
count a spin modulation in space, as has also been stud- 
ied for the one-band Hubbard model 35-40 ' 67 and the t-J 
model. 68-71 We can introduce a stripe in the uniform spin 
density state so that doped holes occupy new levels close 
to the starting Fermi energy keeping the energy loss of 
antiferromagnetic background minimum. The wave func- 
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FIG. 8. Energies of commensurate and incommensurate SDW 
states on 16 x 4 lattice at 8 = 1/8 for Ud = 8. Circles and 
triangles are for 4-lattice and 8-lattice stripes, respectively. 
Squares denote energy for commensurate SDW state. From 
the top e p — 6d — 1.2, 2.0 and 2.4. We impose antiperiodic 
boundary condition in ^-direction and periodic boundary con- 
dition in y-direction. Monte Carlo statistical errors are within 
the size of symbols. 



tion with a stripe can be taken of the Gutzwiller type: 
?P stripe = Pg^ 'stripe- ^stripe is tne Slater determinant 
made from solutions of the Hartree-Fock Hamiltonian 67 
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FIG. 9. Spin (a) and charge (b) densities for incommensurate 
state at 8 = 1/8 for t pp = 0.4, Ud — 8 and e p — td = 2 on 
16 x 4 lattice. The boundary conditions are same as in Fig. 8 



where H® p is the non-interacting part of the Hamiltonian 
H with variational parameter e p — id- (n-di) and (mi) are 
expressed in terms of modulation vectors Q s and Q c for 
spin and charge part, respectively. Including the constant 
part of (ridi) in the definition of variational parameter id, 
we diagonalize the following one-particle Hamiltonian to 
determine ip° stripe : 



H st „pe = H Q dp + ^2\5n di - cr(~l) x * 

ia 

Sridi and m, are assumed to have the form 
Sridi = a/cosh((xi - xf r )/£ c ), 



'mi}d}d ia . 



(16) 



(17) 



(18) 



with parameters a, m, ^ c and t; s 



denotes the posi- 



tion of a stripe. In actual calculations we set £ c — 1 and 
£ s = 1 since the energy expectation values are mostly in- 
dependent of £ c and £ s . Since any eigenfunction of the 
Hamiltonian H® p can be a variational wave function, we 
optimize a instead of fixing it in order to lower the energy 
expectation value further. It is also possible to assume 
that Sridi and rrii oscillate according to cosine curve given 
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FIG. 10. Spin structure function for incommensurate state at 
5 = 1/8 for t pp = 0.4, 0.1 and e p — e d = 2 on 16 x 4 lattice. We 
set Ud = 8. The boundary conditions are same as in Fig. 8. 
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as cos(4ir5xi) and cos(27r5a;.;), respectively. Both meth- 
ods give almost the same results within Monte Carlo sta- 
tistical errors. 

Recent neutron scattering experiments suggested that 
modulation vectors are given by Q s = (n ± 2tt5, it) and 
Q c = (zLAttS, 0) in the underdoped region, where <5 de- 
notes the doping ratio. Here we define n-lattice stripe as 
an incommensurate state with one stripe per n ladders 
for which Q s is given by Q s — (it ± 7r/rt, n). Then the 
incommensurate state predicted by neutron experiments 
for 5 = 1/8 is given by 4-lattice stripe. For the three-band 
model, the transfer t pp between oxygen orbitals plays an 
important role to determine a possible SDW state. If t pp 
is very large, the uniform SDW state is expected to be 
stabilized because holes doped on oxygen sites can move 
around on the lattice producing disorder effect on spin 
ordering uniformly. For small t pp the stripe states are 
considered to be realized. 41 Our motivation to consider 
non-uniform states for the three-band model lies in the 
idea that the distance between stripes may be dependent 
upon t pp , i.e. for small t pp the distance between stripes 
is large, for intermediate values of t pp the 4-lattice stripe 
state is realized and for large t pp the uniform state or 
normal state is stable. 

In Fig. 8 we show the energies for commensurate and 
incommensurate SDW states on 16 x 4 lattice at 6 = 1/8 
as a function of t pp , where we impose the antiperiodic 
and periodic boundary conditions in x- and y-direction, 
respectively, so that the closed shell structure is followed 
for doped holes. We assumed that e p — = 1.2, 2 and 
2.4. The 8-lattice stripe state for small t pp changes into 
uniform state as t pp increases. It shows that incommensu- 
rate states become stable for large level difference e p — td- 
The spin and charge densities of incommensurate state 
are shown in Fig. 9 for t pp = 0.4 and e p — e</ = 2 where 
the charge density is a sum of hole numbers on d-, p x - 
and p^-orbitals at site L. Spin density S z (i) = ndi] — ndi[ 
vanishes at the positions of stripes associated with peaks 
of hole density. The spin structure factor S z (q) really 
has incommensurate peaks as is shown in Fig. 10. The 
Figures 11(a) and 11(b) present the energies of incom- 
mensurate states for 16 x 16 lattice (which contains 768 
atoms) where we set antiperiodic and periodic boundary 
conditions in x- and y-direction, respectively, for (a) and 
in y- and x- direction, respectively, for (b). Both figures 
give almost the same results as an evidence that the ef- 
fect of boundary conditions is small for 16 x 16 system. 
As expected, the structure of incommensurate state is 
dependent upon the values of t pp . 

Let us turn to a discussion of the energy gain due to 
a formation of stripes, which is estimated from an ex- 
trapolation to the bulk limit as shown in Fig. 12. One 
notes that the energy gain increases as the system size 
increases. The energy gain per site for 4-lattice stripe 
state is given by ~ 0.015id p ~ 22.5meV. Furthermore the 
energy difference between commensurate and incommen- 
surate states is found to be finite in the bulk limit, which 
is shown in Fig. 13. Thus within VMC the stripe state 
with spin modulation is stable at 8 = 1/8 doping. 
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FIG. 11. Energies at 5 = 1/8 for 16 x 16 lattice. Parameters 
are given by e p — = 2, Ud = 8 and t pp — 0.4. Symbols are the 
same as in Fig. 8. For (a) boundary conditions are antiperiodic 
and periodic in x- and {/-direction, respectively, and for (b) 
periodic and antiperiodic boundary conditions are imposed in 
x- and {/-direction, respectively. Monte Carlo statistical errors 
are within the size of symbols. 



The antiferromagnetic order parameter m in eq.(18) 
is of the order of 0.5^ ~ 0.75eV, while the SC order 
parameter A s (which gives the minimum of energy) is 
of the order of 0.01 - 0.015^ = 15mcV- 20meV at 
S ~ 0.2. The magnitude of SC order parameter agrees 
with measurements of tunneling spectroscopy 72 ' 73 where 
A s is estimated as A s ~ 17meV for YBCO sample. 72 
The antiferromagnetic order parameter is larger than SC 
order parameter at least by one order of magnitude. The 
charge order parameter a in eq.(17) is small and negligi- 
ble compared to mj. 



V. Summary 

We have presented our evaluations for the 2D three- 
band Hubbard model based on the variational Monte 
Carlo method. Our work is regarded as a starting step 
for more sophisticated calculations in future such as the 
inclusion of correlation factors of Jastrow type or Green 
function Monte Carlo approaches. The SC energy scales 
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FIG. 12. E norma i — E incom as a function of l/N. Parameters 
are given by t pp = 0.4, Ud = 8 and e p — = 2. Solid circles 
are for 24 x 6, 32 x 8 and 40 x 10. Open circle is for 16 x 16. 
We set antiperiodic and periodic boundary conditions in x- 
and y- direction, respectively. 
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FIG. 13. E com — Eincom as a function of l/N. Parameters are 
given by t pp = 0.4, Ud = 8 and e p — = 2. Circles are for 
rectangular lattices and the square is for 16 x 16. Boundary 
conditions are the same as in Fig. 12 



obtained from our evaluations are consistent with exper- 
imental indications, which provides a support to our ap- 



According to VMC the attractive interaction works 
for d-wave pairing due to electron correlations. The 
strength of Ud is also important to determine the phase 
boundary of the SDW phase. If Ud is extremely large, 
the SDW region extends up to large doping for which 
the d-wave region is restricted to infinitesimally small 
region near the boundary of antiferromagnetic phase. 
For intermediate values of Ud and e p — td the SDW re- 
gion is reduced and the d-wave superconducting phase 
may exist. The fact that the SC condensation energy 
agrees reasonably with the experimental data for opti- 
mally doped samples supports our computations. The 
magnitude of SC order parameter is also consistent with 
tunneling spectroscopy experiments. From our data for 
AE.sc and A s and the relation N(0)A 2 J2 = AE S c, 
the effective density of state N(0) can be estimated as 
AT(0) ~ 3 ~ 6.7(eV)~ 1 ~ 4.4 - 10/t dp at S ~ 0.2 in the 
overdoped region, which is not far from the BCS estimate 
2V(0) - 2 to S(eV)- 1 by using Y(0)(fc B T c ) 2 /2 for opti- 
mally doped YBCO. 64 We expect that the pure d-wave 
state from optimal to overdoped regions may be described 
by the projected-BCS wave function. The phase diagram 
for electron-doping is consistent with the available ex- 
perimental indications suggesting that the properties of 
electron-doped materials may be understood within our 
approach. In the SDW region the incommensurate spin 
structures are stabilized for the low-doping case to keep 
the energy loss minimum due to disorder effect caused 
by holes. A competition among the uniform SDW state, 
SDW state with stripes, and pure d-wave SC is highly 
non-trivial. A picture for the hole-doping case followed 
from our evaluations is that a stripe state is stable in the 
underdoped region and changes into the d-wave SC in the 
overdoped region. 
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